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We consider a local quench where two free-fermion half-chains are coupled via a defect. 
We show that the logarithmic increase of the entanglement entropy is governed by the same 
effective central charge which appears in the ground-state properties and which is known 
exactly. For unequal initial filling of the half-chains, we determine the linear increase of the 
entanglement entropy. 

I. INTRODUCTION 

The entanglement across defects in critical quantum chains is interesting, because in free-fermion 
systems it varies continuously with the defect strength For two half-chains of length L, one 
finds a logarithmic behaviour of the von Neumann entanglement entropy 

S = ^lnL (1) 

where < c c e < c with c (equal to 1/2 or 1) denoting the central charge of the model. Since e s , 
which is the effective number of terms in the Schmidt decomposition, then follows a power law, 
one can view c c q as a variable critical exponent. It was studied in a number of papers for discrete 



nri 



n 



pHa] and continuous [6H8|] systems and an exact analytical formula was obtained 3j, also for the 
Renyi entropy 0, [ill and for bosons [{], • Technically, the variation is caused by a gap in the 
single-particle eigenvalue spectrum of the reduced density matrix (RDM) and the relevant physical 
parameter is the transmission amplitude of the defect. 

A similar situation is found if one considers a local quench where one connects two initially 
separated half-chains by a defect. The entanglement entropy then increases in time as 

S= d -flnt (2) 



as long as t <C L for finite L and for all times if L is infinite. This was found either 



numerically Q| or by evaluating expressions based on (lie counting statistics numerically 



rjurely 

H 
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No explicit formulae for c e g have been given so far. However, it was observed in the numerics of 
the transverse Ising model that c e ff equals c e g, as in the the homogeneous case, where the same c 
appears in ([I]) and ([2]) and the formulae follow from conformal considerations 141 ] . 

In this note, we look at this problem once again and study the evolution after the quench for an 
XX (hopping) model via the time-dependent single-particle RDM eigenvalues. We show that those 
for the defect case are related by a simple exact formula to those for the homogeneous case, if the 
defect is conformal (i.e. scale- free). The relation is the same as in the static case and thus implies 
again a gap in the spectrum. This allows to take over the ground-state results and thereby proves 
the equality of c e g and c e ff. Since c c ff is exactly known, it also provides an analytical formula for 
c e ff. One can check that it fully agrees with the expressions and results of the counting statistics. 

We also consider the case where the initial filling of the two half-chains is different. Then a 
somewhat different relation for the spectra exists where the effect of backscattering is manifest. 
The current which sets in after the quench leads to a steady generation of entanglement between 
transmitted and reflected parts of the wavefunction and thus to a linear increase of the entropy with 
time. We give analytical expressions for the coefficient both for conformal and for non-conformal 
defects. 



II. MODEL AND METHOD 



We study free fermions hopping on a finite open chain of 2L sites. The time evolution for t > 
is governed by the Hamiltonian 

1 2L 

H = ^ ^2 H-m,n c m c n (3) 
m,n=l 

where the nonzero matrix elements are 

/ / J-l m±L , , , 

I —A m = L 

In the middle of the chain there is a defect in the form of a modified bond supplemented with site 
energies on both sides and characterized by the parameter A. It will be referred to as the conformal 
defect and is used to derive all the exact relations. In addition, we will also consider simple bond 
defects without site energies. The initial system is composed of two disconnected half-chains and 
the initial state will be specified in the corresponding sections. 

In the homogeneous case, A = 1, the Hamiltonian (J3j) is diagonalized by a Fourier transform 
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and the eigenvectors and eigenvalues are given by 



a ( \ / 2 • ^ km n nk fK\ 
<Pfc(?7i) = A/— sin— , ilh = — cos— 5) 

YK ' V2L + 1 2L + 1' 2L + 1 KJ 

where = 1, . . . , 2L. For the conformal defect the solutions are related to those of the homogeneous 
chain via 

, I OLk<pk{fn) 1 < m < L f 

4> k (m) = < , n k = n k (6) 

I /3fc0fc(w<) L < m < 2L 
and thus the eigenvectors are simply rescaled on the left and right hand side of the defect with the 
scaling factors 



l + (-l)Vl-A 2 , ft 2 = l-(-l)Vl-A 2 (7) 



Note that, apart from the alternation, a k and j3 k are independent of k and lead to a constant 
transmission coefficient T = A 2 . Therefore one has here a lattice realization of the special scale- 
free defects which can be constructed by gluing together two conformal field theories [lol . [lfj . The 
bond defect is non-conformal and leads to a transmission which depends on the wavelength of the 
incoming particle. 

For the entanglement between left and right halves, one needs the RDM for a half-chain which 
has the form p = e~^/Z with a free-fermion effective Hamiltonian H, see [if]]. In the following, 
the single-particle eigenvalues of T~L will be called 2u>;(i). They are related via 

4® = ^WTT (8) 

to the eigenvalues of the correlation matrix C'(t) which is the full matrix C'(t) with elements 
{cln{t)c n (t)) restricted to 1 < m,n < L. For A = 1, the quench will be called homogeneous and the 
corresponding quantities are denoted by 2ei (t) , Q (t) and C (t) , respectively. 

Therefore one has to find the time-dependent correlation matrix which can be written in the 
Heisenberg form 

C'(i) = e ifl '*C(0)e- ifl '* (9) 

where H' denotes the matrix with elements H' mn and C(0) contains the initial correlations. The 
entanglement entropy is then given by 

S{t) = £ ln(l + e-**«>) + £ J^- = £ H(t>(t)) (10) 
I i i 

where H{x) = — xlnx — (1 — x) ln(l — x). The first (second) form proves to be useful for quenches 
from equal (unequal) fillings. 
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III. QUENCH FROM EQUAL FILLINGS 

In this case, the initial correlation matrix has the block form 

, C° , 

C(0) = (11) 

c° 

where the L x L matrix C° refers to a half-chain in the ground state and has elements 

L 

C^ n = J24H4(n)n k (12) 

k=l 

Here n k is the occupation number of mode k and the wavefunctions 4^.(m) have the form (0) with 
2L replaced by L. To evaluate ([9]) one expresses the exponential operators as 

2L 

L±m't\ =£#( TO )#(n)eF«n** (13) 

\ / van 

k=l 

Using @, defining the overlap matrix 

L 2L 

B kl = Y J Mj)$ti) = (-l) (k+l) E M)<fi(j-L) (14) 

3=1 j=L+l 

and noting that a k a k i + (— l) fc+fc PkPk' = 2 for k — k' even and zero otherwise, one arrives at 

2L L 

C' mn {t)= 2 E Yj «fc«fc'^fc"^'fc"^(^)0fc'(n)^"e i( ^~ nfc,) * (15) 
k,k'=l k"=i 

k—k'even 

The only dependence on the defect strength A is in the factors a k ay . The point now is that one can 
find a simple connection with the homogeneous quench if one considers the matrix (2C'(i) — l) 2 . 
Then a further overlap matrix appears with elements 

L I 



2(2L + 1 



sinf(Jfe-i) sinf(fc + Z) 

7r(fc-Q i 7r(fc+Q 

S111 2(2L+1) Sm 2(2L+1) 



(16) 



which, apart from the diagonal ones, vanish for k — I even. The matrix elements of (2C'(t) — l) 2 
are 6-fold sums where each summand is proportional to a k a k iaia\i A k 'i with k — k' and / — I' even. 
For k' — I odd, this gives a k a k iaiai' = A 2 and thus a simple prefactor. For k' = I, one has An = 1/2 
and a summand 

a k a\a v B kk nB lk nB ll »B in ^k(ra)(\) V {n)(2n k » - l){2n v , - l)^ *""'')' (17) 
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where a k a^ay = 2 — A 2 + (—l) l 2y/l — A 2 for the allowed indices. The remaining sums can then be 
evaluated using the identities 

2L 2L L 

y BiyiBu'i = 5 k "i", Bik" Bu» = 0, B kk nBi' k n = A kV (18) 

1=1 i=i fc"=i 

Summing over I therefore eliminates the alternating terms and enforces k" = I" which leads to 

(2n k /> — l) 2 = 1 independent of k" and of the filling. The final sum over k gives (1 — A 2 /2)5 mn . In 

this way one arrives at the formula 



(2C'(t) - 1)L = A 2 (2C(t) - 1)L + (1 - X 2 )5 n 



(19) 



As a consequence, the single-particle eigenvalues of the corresponding RDMs are related by 
tanh 2 QJi(t) = A 2 tanh 2 £i{t) + 1 — A 2 , or alternatively, writing A = s for the transmission amplitude, 
by 



ch^(t) = -ch Ei(t) 



(20) 

This is ex-nelly I he .same relation one finds lor the sl-al-ic defect problem ■>] and gives a gap in the 
a;;-spectrum. Note that (|19p can also be written as C'(t)(l — C(t)) = A 2 C(t)(l — C(t)) and is then 
identical to the relation for the overlap matrix A in a (static) continuum system [6H8J. 
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FIG. 1: Low- lying eigenvalues 2ei(t) after homogeneous local quench with 2L = 400 (points) and 2ei(£) of 
an interval in an infinite system (lines). The data is shown for even I and for half- filling. 



The problem is now reduced to that of the homogeneous quench, but one still needs the £i(t). 
The known results give a logarithmic behaviour of S(t) with a prefactor 1 /3 that is similar to the 
scaling of the equilibrium entropy of an interval t in an infinite chain 11411 . In the static case, one 
knows that the lowdying £i(£) have a spacing 7r 2 /21n^ for large ln£ [171 |. This gives a density of 



6 



states proportional to \n£ and the logarithmic variation of S(i). The similarity of the conformal 
results for the entropy suggests that the £i(t) have an analogous behaviour. This is in fact the 
case and is shown in Fig. [TJ Plotted are the lowest £i(t) in a homogeneous quench together with 
the lowest £[(£) for a segment of length I in an infinite chain. Apart from some oscillations, they 
coincide closely. 

With this observation, one can now follow the same steps as in the static case and write S in 
(fTTJj) as an integral over e. In this way one finds c e g = c e g, where c e g is 12/ it 2 times the integral 
J(s) given in (26) of [3j. It is depicted on Fig. EJboth as a function of s and T = s 2 , to allow 
simple comparison with the numerical results of [111 ]. The function is nonanalytic around s = 0. 
For the simple bond defect, one knows from the static case that the parameter s is the transmission 
amplitude at the Fermi level. The relation between oj and e is then only valid as the levels become 
dense. 
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FIG. 2: c e fF as a function of s (solid line) as well as T — s 2 (dashed line). 



IV. QUENCH FROM UNEQUAL FILLINGS 

If the initial fillings are unequal, a steady particle-current results after connecting the half- 



chains. In the homogeneous case, the time evolution of the density profile 181, the particle number 

□ □ 
fluctuations in a half-chain [19] as well as the entanglement entropy [20] have been studied previ- 
ously. The fluctuations and the entropy both grow logarithmically in time. 

We first consider the case where the left hand side of the chain is completely filled and the right 
hand side is empty. The initial correlation matrix then reads 

C(0): | 1 ° | (21! 
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where 1 denotes the L x L identity matrix. Expanding into eigenvectors one has 

2L 

C' mn {t) = ola 2 ^ k ( m )A kl Un)e mk - ni)t = X 2 C mn {t) + (1 - X 2 )5 n 



(22) 



k,l=l 



where in the second step we used the property a^af = A for all A: — Z odd, where the overlap 
matrix A is nonvanishing. For k = I, the alternating piece in a\ gives zero upon summation and 
the constant piece leads to the second term on the right hand side. Remarkably, one finds the same 
relation as in Eq. (|19p but now for the correlation matrix itself and thus also for its eigenvalues 



c;(o = A 2 ow + i-A 2 



(23) 



The eigenvalues Cz(^) and Q(t) are shown on the left of Fig. [3] and have a simple interpretation. 
As the half-chains are connected, particles leave the left part of the system and a steady current 
results. In the pure case, this implies a growing number of Q(t) = eigenvalues. In the presence 
of a defect, however, there is a probability R = 1 — T for backscattering and the corresponding 
eigenvalues are C;(^) = 1 — A 2 = i2. It is interesting to compare with the case of a simple bond 
defect, which is done on the right of Fig. [3J Here such a simple relation as (|23p does not hold, 
since the transmission coefficient depends on the momentum q of the particle as 



sin 2 q 



ch v — cos 2 q 



(24) 



where we defined A = e v and < q < it. Because of this variation, the C/W curves are not flat. 
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FIG. 3: Eigenvalues Q(t) (open symbols) and (filled symbols) for different times and 2L = 200 sites. 
Left: conformal defect. Right: bond defect. 



The common feature of both defects is that the number of nonzero Cz(^) eigenvalues, which have 
a finite contribution to the entropy, grows with time proportionally to the number of transmitted 
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particles. Because of the steady flux in the center of the chain, this results in a linear contribution 
to the entanglement. Numerically, this can be observed even for coupling strengths arbitrarily 
close to A = 1, where the entropy becomes logarithmic, S ~ 1/6 log t 
well fitted with the ansatz 



201 ] . In general, S can be 



S(t) =crf + /?lnt + 7 



(25) 



and the coefficients are shown in Fig. 2] In the conformal case, the linear part can be calculated, 
using the second form of the entropy in Eq. ()10p . as H(X 2 ) t/ir. This is just the contribution of a 
single £>(i) on the flat part of the spectrum multiplied by the length of the plateau, which equals 



181 ] . This result, shown by 



the total number t/ir of transmitted particles in the homogeneous case 
the solid line on the left of the Fig. H] agrees perfectly with the fitted values of a. Furthermore, it 
also agrees with the results of 



11 



13] for the entropy evolution in a quantum point contact in the 



high-bias regime. 
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FIG. 4: Fitted coefficients of the ansatz S(t) = at + /31n£ + 7 together with the analytical formulas obtained 
for the slope a. Left: conformal defect. Right: bond defect. 



In the non-conformal case, the linear growth of the entropy can be obtained by distinguishing the 
different momentum states and associating a factor H(T q ) with the entanglement which is produced 
between transmitted and reflected parts of the wavefunction. Thus one writes for L — > oo 



a 



(26) 



where v q = dVt q /dq = sing is the velocity of the incoming particle and measures the flux. For the 
conformal defect where T = A 2 one obtains the previous result. For the bond defect, the integral 
can be carried out by substituting Eq. (|24p and using T q as integration variable. A lengthy 
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calculation then yields a = I(s)/tt with 
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(27) 



4s 2 j s ~ \ 1 + s / 

where Lia(s) denotes the dilogarithm function and s = 1/chi^ = 2/(A + l/A) is the same amplitude 
that appears in the equilibrium formulae. This analytical result is again in perfect agreement with 
the results of the data fits, as shown on the right of Fig. HJ 

Note, that the structure of Eq. (|26l) is very similar to the one found by Fagotti and Calabrese 
for a global quench in the XY model 2l| . However, in the present case entanglement is created 
only locally (but steadily) at the defect, in contrast with the global quench where entangled quasi- 
particle pairs are created everywhere but only at t = [22]]. It would be interesting to check if this 
semiclassical picture can be used to calculate other relevant observables as was found recently for 
the global quench Q] • 

The above arguments can be generalized to arbitrary initial filling factors ni and n r on the 
left and right. Assuming ni > n r , the single-particle states with q < n r Tr will be filled and thus 
balanced on both sides. In terms of the the number of eigenvalues on the flat part of the 

spectrum decreases and one has a more general relation which interpolates between (|23p and (|20p . 
In the semiclassical picture only wavenumbers n r Tr < q < niir will contribute to the current and 
the integral in Eq. (|26p has to be carried out only on this interval. As shown in Fig. [5l this gives 
again very good agreement with the numerical data, further supporting the semiclassical picture 
described above. It is possible to evaluate the integral analytically also in this case. 
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FIG. 5: Slope a of the linear part of the entropy for several filling factors as obtained by fits to the data 
(points) as well as evaluating the integrals (lines) 
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V. CONCLUDING REMARKS 

There is a close connection of our findings with the work on counting statistics by Klich and 
Levitov 111 ]. For the quench with equal fillings, they obtain a density of states for the eigenvalues 



£' which vanishes in an interval around £' = 1/2. Thus the spectrum has a gap, and this is just the 
situation described by (|20p . Moreover, if one converts their entropy expression (10) into an integral 
over e one finds exactly the integral for c e fr evaluated in [3|. Finally, the relation between their 
parameters A and A* is another form of the dispersion relation (|2Uj) derived here. Our approach, 
focussing on the eigenvalues themselves, is more direct and shows the connection with the static 
defect problem very clearly. The only open point is an analytical derivation of the £/(i)-behaviour, 
which we inferred from confer mal results and numerics. 

In this context one should mention t hat we discussed only the case t -C L. The general CFT 
formula for the homogeneous quench is 



25] 



S(t) = °- In 



2L TTVFt 

sin 



7T 2L 



+ const. (28) 



where the Fermi- velocity vp is set by the filling of the half-chains. This results again from the 
scaling of the £i(t) with the logarithmic factor and the formula for the defect case is obtained by 
substituting c — > c e ff. 

The relation c e fr = c e fr also holds for the corresponding coefficients in the Renyi entropies S n 
which are rather simple for 52 and 53 . Similarly, the considerations for the biased quench 

can be adapted to the Renyi functions. 

Finally, a quench from unequal fillings also arises in continuum models if one removes a wall 
which initially confines the particles. A detailed study of the entanglement evolution has appeared 
recently 



261 ] and can probably be extended to the defect case. 
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